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ABSTRACT 

We present the results from a Chcmdra pilot study of 12 massive mergers selected from Galaxy Zoo. 

The sample includes major mergers down to a host galaxy mass of 10 11 Mq that already have optical 
AGN signatures in at least one of the progenitors. We find that the coincidences of optically selected 
active nuclei with mildly obscured (iV// < 1.1 x 10 22 cm -2 ) X-ray nuclei are relatively common (8/12), 
but the detections are too faint (< 40 counts per nucleus; / 2-10 kev < 1-2 x 10“ 13 erg s -1 cm -2 ) to 
separate starburst and nuclear activity as the origin of the X-ray emission. Only one merger is found 
to have confirmed binary X-ray nuclei, though the X-ray emission from its southern nucleus could be 
due solely to star formation. Thus, the occurrences of binary AGN in these mergers are rare (0-8%), 
unless most merger-induced active nuclei are very heavily obscured or Compton thick. 

Subject headings: galaxies: active — X-rays: galaxies 


1. INTRODUCTION 

Major mergers are a key component of current mod- 
els for galaxy formation in a ACDM Universe. Merg- 
ers can disrupt the star-forming gas and stellar disks of 
the progenitors, trigger a powerful burst of star forma- 
tion, and reshape the remaining stellar content into a 
bulge. Perhaps with a small time delay, the supermas- 
sive black holes may feed on gas from the destabilized or 
destroyed disk, injecting energy in the form of radiation 
or kinetic outflows that sweep the remnant clear of dust 
and gas. First proposed by Sanders et al. (1988), this 
picture directly finks the triggering of active galactic nu- 
clei (AGNs) phases to the destructive potential induced 
by a merger. Recent semi-analytic models and hydrody- 
namic simulations have adopted this scenario to explain 
the fueling of AGNs and the red spheroidal remnants that 
are difficult to reproduce without some kind of “AGN 
feedback” (Springel et al. 2005; Di Matteo et al. 2005; 
Hopkins et al. 2006, 2008; Somerville et al. 2008). 

1 Observational Cosmology Laboratory, NASA/GSFC, Green- 
belt, MD 20771, USA; stacy.h.teng@nasa.gov 

2 Department of Astronomy, University of Maryland, College 
Park, MD 20742, USA 

3 Department of Physics, Yale University, New Haven, CT 
06511. USA 

4 Yale Center for Astronomy and Astrophysics, Yale Univer- 
sity, PO. Box 208121, New Haven, CT 06520, USA 

® Department of Astronomy, Yale University, New Haven, CT 
06511, USA 

8 Department of Physics, University of Oxford, Keble Road, 
Oxford, OX1 3RH, UK 

7 Department of Astronomy, Yonsei University, Seoul 120-749, 
Republic of Korea 

” The Harriet W. Sheridan Center for Teaching and Learning, 
Brown University, P.O. Box 1912, Providence, RI, 02912, USA 

9 Department of Physics & Astronomy, 206 Gallalee Hall, 514 
University Blvd., University of Alabama, Tuscaloosa, AL 35487- 
034, USA 

19 Departamento de Astronomia, Universidad de Concepcion, 
Casilla 160-C, Concepcidn, Chile 

11 NASA Postdoctoral Program Fellow 

12 Einstein Fellow 


In principle, major mergers carry with them two black 
holes, both of which may be accreting and be visible 
as distinct AGN during a phase of abundant gas avail- 
ability that a major, gas-rich merger represents. Yet the 
evidence associating AGN phases with major mergers re- 
mains contested (De Robertis et al. 1998; Malkan et al. 
1998; Schmitt 2001; Pierce et al. 2007; Georgakakis et al. 
2009; Gabor et al. 2009; Schawinski ct al. 2011). Large 
optical surveys using Sloan Digital Sky Survey (SDSS) 
data have found ~3.6% of spectroscopically confirmed 
AGNs are in closed binaries (~5-100 kpc separation; Liu 
et al. 2011). The DEEP2 survey also found that binary 
AGN exist in ~2.2% (2/91) of red galaxies with type 
2 Seyfert optical spectra (Gerke et al. 2007; Coraerford 
et al. 2009). However, optical surveys can easily miss 
obscured AGNs especially in merger systems where the 
gas is driven toward the center through dissipation (e.g., 
Hopkins et al. 2008). X-ray surveys are needed to identify 
the more highly obscured systems (N fI > 10 20 cm -2 ). 
We know of only a small number of binary AGN resolved 
directly using X-ray observations (e.g., Komossa et al. 
2003; Guainazzi et al. 2Q05; Hudson et al. 2006; Bianchi 
et al. 2008; Foreman et al. 2009; Comerford et al. 2011; 
Fabbiano et al. 2011). The intrinsic frequency of binary 
AGN phases has not been observationally constrained, as 
the separation of individual X-ray sources is not possible 
in high-redshift sources and there has been no systematic 
search for such systems in known mergers. A study of 
the hast galaxies of 185 nearby (z < 0.05) BAT AGNs by 
Koss et Al. (2011) found that these X-ray selected AGNs 
are preferentially found in massive galaxies with large 
bulge-to-disk ratios and large supermassive black holes. 
This may imply that the frequency of binary AGNs are 
higher in massive mergers. 

In order to quantify the intrinsic frequency of double 
AGNs in the local universe, we embarked on a study of 
the presence of binary AGNs and their dependence on the 
mass of the host galaxies. The results from the present 
survey represent a pilot effort as the sample is comprised 
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of only the most massive galaxies in the Galaxy Zoo 
merger sample. Thus, this paper aims to quantify the 
intrinsic frequency of double AGNs in the mass limit 
down to ~ 10 n Mq using a study of 12 merging galaxies 
with the Chandra X-ray observatory. Throughout this 
paper, we adopt H o = 71 km s -1 Mpc -1 , Ojf = 0.3, and 
fl A = 9.7. 

2. SAMPLE SELECTION 

The initial parent sample was created from 3003 merg- 
ers identified via visual inspection by citizen scientists 
taking part in the Galaxy Zoo project (Lintott et al. 2008; 
Darg et al. 2010a, b). To date, this is the largest unbi- 
ased sample of mergers in the local Universe (0.005 < 
z < 0.1). From this catalog, we whittled the sample 
down to only 328 major mergers (i.e., with a mass ra- 
tio of 3:1 or less) in which the primary galaxy has a 
mass 13 > 10 n Mg; we also required that the galaxies 
have SDSS spectroscopic data with signal-to-noise ratio 
greater than 3 and that at least one of the nuclei shows 
an AGN signature based on [O IIIJ/H/J and [N II] /Ha 
narrow line ratios (Baldwin et al. 1981; Kewley et al. 
2006). In addition, the AGN must be relatively lumi- 
nous, with L[ 0 in] > 10 41 erg s _1 (yellow points in Fig- 
ure 1) This last cut is to ensure that the AGN will 
be luminous enough for detection in the X-rays. Most 
of the objects in this final sample are LINERs or AGN 
and star-forming composites, from which we selected 12 
that spanned the full starburst-composite-LINER/AGN 
range. The primary nuclei of the selected sample (yellow 
dots circled in red) also have roughly even distributions 
in the composite but not extreme starbursts (40%) and 
the LINER (53%) regimes of Figure 1 and these mergers 
have projected nuclear separations between 5 to 14 kpc. 
The optical line ratios are taken from Oh et al. (2011). 
The 12 galaxies in our sample are listed in Table 1; 
for convenience, we will identify the targets as Galaxy 
Zoo (GZ) objects throughout this paper. For clarity, we 
will refer to the merging galaxies as either “mergers” or 
“galaxies” and the individual progenitors of these merg- 
ers as “nuclei” for the remainder of this paper. Figure 2 
is a collage of the SDSS images of the sample. All of 
these mergers show disturbed morphology indicative of 
tidal disruptions. 

3. OBSERVATIONS, DATA REDUCTION AND ANALYSIS 

The 12 mergers were observed with Chandra between 
November 2010 find May 2011 (PI: Schawinski). Each 
pair of galaxies was observed in a single exposure of 
4.9 ks. For each merger, the more massive primary nu- 
cleus was placed at the nominal aim point of the ACIS-I3 
chip.The progenitors of the mergers are close enough that 
the secondary nuclei were also within the same chip. 

The data were reduced using CIAO version 4.3 and 
CALDB version 4.4.3. The data reduction followed pro- 
cedures outlined in the Science Analysis Threads for 
ACIS imaging data on the CIAO web page 14 . For each 

13 The stellar masses of the galaxies were calculated following 
the methodology outlined in Schawinski et al. (2010). Briefly, mea- 
surements from the five SDSS photometric bands were fitted to a 
library of model star formation histories generated from Maraston 
(1998, 2005) stellar models. Stellar masses are measured by finding 
the minimum of the x 2 statistic in the parameter space probed. 

14 http:// cxc.harvard.edu/ciao/threads/index.html 


of the 12 detected nuclei, we extracted counts in the soft 
(0.5-2 keV) and hard (2-8 keV) bands. Two of these 
nuclei were detected having only two counts in the soft 
band, but none in the hard band. Therefore, only 10 
nuclei have valid hardness ratios (HRs 15 ; Table 1) for es- 
timating spectral properties. For the rest of this paper, 
only the 10 nuclei with valid HRs will be discussed. 

An HR analysis was performed using the latest version 
of FTOOLS released as part of HEASoft 6.11. Due to 
the low number of counts measured from our sample, we 
followed Teng et al. (2005) and used the measured HR 
to estimate a photon index (T) by assuming a redshifted 
power law model modified only by Galactic absorption. 
The nominal HR and the estimated photon indices are 
tabulated in Table 1. The X-ray fluxes were then esti- 
mated using PIMMS by assuming the HR-derived photon 
indices and the count rates from the observations. 

4. DISCUSSION 

The shapes of the X-ray spectra differ for AGNs, ob- 
scured AGNs, starbursts, and AGN-star forming com- 
posites. Typically, unobscured AGNs have spectra that 
are well-represented by a power law with photon index 
of ~1.8. Obscurations affect the lower energy (<2 keV) 
photons more readily than the higher energy photons and 
thus flatten or harden the AGN spectra. Starburst spec- 
tra are dominated by emission in the lower energies, but 
low-mass X-ray binaries tend to have relatively flat spec- 
tra. Composite objects tend to have softened spectra 
compared to simple AGN spectra. 

Given that these mergers contain optically selected 
AGNs, it is unsurprising that eight of the 10 nuclei have 
HRs that are consistent with the canonical spectral shape 
of unobscured AGNs (r ~ 1.7 — 2.1). As many as five 
could be steeper (GZ IE, GZ 4S, GZ 5S, GZ 7SW, and 
GZ 9S), as if star formation is a significant contribu- 
tor, though the errors in HR and T allow for unobscured 
AGN values. Similarly, three nuclei (GZ 2N, GZ 9N, 
and GZ 11S) have nominally flat spectra, implying dom- 
inance from star formation or obscured nuclear activity. 
However, the errors in HR cannot rule out unobscured 
AGN as the source of the X-ray emission. Finally, two 
nuclei (GZ 3S and GZ 10S) have flat or inverted photon 
indices (F < 1.45 after accounting for the measurement 
errors), suggesting some level of obscuration. If we as- 
sume a power law with F fixed at 1.8, the observed HRs 
imply column densities (Nr ) < 10 21-22 cm -2 (Table 1). 
At these column densities, the HR estimates of the 2- 
10 keV luminosity are reliable to within ~40% (Teng 
& Veilleux 2010). These columns do not suggest the 
presence of Compton-thick nuclei though there remains 
a possibility of leaky, Compton-thick absorbers. 

4.1. Starburst Contamination 

Of the 12 mergers in the sample, one has no X-ray 
detection (GZ 8) and only one (GZ 9) exhibits binary 
X-ray nuclei (Figure 3). The remaining 10 mergers have 
one detected nucleus each. In GZ 5, GZ 7, and GZ 12, 
the X-ray-detected nucleus is not the one with an op- 
tical AGN classification, so in that sense they are dou- 
ble nuclei. In addition, the detected southern nucleus 

18 HR = , where H and S are the total counts in the hard 

and soft bands, respectively. 
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of GZ 10 has extended soft X-ray emission (Figure 3), 
suggesting contribution from star formation. This raises 
the question of whether more of the detections might be 
contaminated by star formation. 

To explore this possibility, we compare star formation 
rates derived from the SDSS u band luminosities follow- 
ing Hopkins et al. (2003) with those derived from the 2- 
10 keV luminosity following Ranalli et al. (2003) in Fig- 
ure 4. When compared with the SDSS u-band derived 
star formation rates (Figure 4), the X-ray derived star 
formation rates of four nuclei (GZ IE, GZ 5S, GZ 7SW, 
and GZ 9S) have unconstrained lower limits. While the 
nominal X-ray derived star formation rates are above the 
line of equality implying the presence of AGNs, we can- 
not rule out the possibility that the X-ray emission can 
be accounted for solely by star formation in these four nu- 
clei. The X-ray luminosities of the remaining nuclei are 
above those expected from star formation even after the 
consideration of the <40% error in the calculation of the 
X-ray luminosity, consistent with additional contribution 
to the X-ray luminosity by nuclear activity. Accounting 
for the error bars, the southern nucleus in GZ 9 may also 
be dominated by star formation, suggesting GZ 9 does 
not contain an AGN pair. 

4.2. Compton-thick Nuclei 

Three nuclei of the 12 SDSS-selected mergers are not 
detected in X-rays (GZ 5N, GZ 7NE, and GZ 12E). We 
already know these have optically identified AGN com- 
ponents, so it is unclear whether the non-detections are 
due to faint AGNs (two of the three have the highest 
redshifts in our sample) or Compton-thick AGNs. If 
we assume these are faint AGNs, a power law model 
with T = 1.8 and mild absorption from the Milky Way 
places upper limits to the luminosity of these objects. 
Not accounting for intrinsic absorption, the 2-10 keV lu- 
minosity for GZ 5N is < 5 x 10 39 erg s _1 cm -2 and 
< 5 x 10 4 ° erg s -1 cm -2 for GZ 7NE and GZ 12E. In 
the case of the Compton-thick AGNs, the optical signa- 
ture is coming from the much larger scale narrow- and 
broad-line regions while the X-ray is sensitive to the small 
scale emission from the black hole itself. The presence 
of undetected obscured nuclei would affect our statistics 
of the frequency of binary AGNs. It is unlikely that 
all of the secondary nuclei contain Compton-thick X-ray 
sources, unless an obscured phase, is common to mergers 
(unlike isolated AGNs). Even without a merger-induced 
obscured phase, the number of heavily obscured AGNs is 
comparable to the number of less obscured AGNs (Treis- 
ter et al. 2009); the presence of Compton-thick nuclei 
remains a possibility. 

While the individual detected nuclei have too few 
counts for spectral fitting to definitively establish 
whether Compton-thick AGNs are present, we consid- 
ered the cumulative rest-frame photon distribution of the 
detected nuclei in the hard band. We compared this ob- 
served distribution with the expected photon distribu- 
tions from unobscured AGNs and from Compton-thick 
AGNs. In the former case, we assumed a single unab- 
sorbed power law; in the latter case, we assumed a power 
law with an iron emission line at 6.4 keV with an equiv- 
alent width of 1 keV, a typical signature of Compton- 
thick AGNs. For both cases, the total photon counts 
were normalized to be the same as the total detected 


counts. In Figure 5 we plot the cumulative distribution 
of the detected photons in our sample. There is no clear 
distinction between the observed distribution with either 
model. In fact, the two-tailed Kolmogorov-Smirnov (K- 
S) test statistics for the two cases are nearly identical. 
As a sanity check, we compared the two modeled distri- 
butions with each other and there is a clear difference 
at the ~80% confidence level. Therefore, we cannot rule 
out the possibility that Compton-thick AGNs are present 
at the level that we are able to detected these sources. 

4.3. JJ/LIRGs in Formation? 

In theoretical models of galaxy mergers (e.g., Hopkins 
et al. 2008), luminous and ultraluminous infrared galax- 
ies (U/LIRGs) represent a stage that mergers go through 
before the formation of elliptical galaxies. Initially, tidal 
torques enhance star formation and black hole accretion. 
Then in the final coalescence of the galaxies, massive in- 
flows of gas trigger starbursts with strengths similar to 
those inferred for U/LIRGs. 

The mergers in our sample appear to be the predeces- 
sors to U/LIRGs in this evolutionary picture. The X-ray 
luminosities estimated for our mergers are approximately 
10 times lower than those observed in most U/LIRGs, 
but are consistent with the lower end of the range mea- 
sured in LIRGs (Teng & Veilleux 2010; Lehmer et al. 
2010; Iwasawa et al. 2011). This implies mergers in our 
sample are in the earliest stages of interaction, where the 
growth of the central black hole has not yet peaked. 

The incidence of binary AGNs in U/LIRGs is also rare. 
The Revised Bright Galaxy Survey (RBGS; Sanders et 
al. 2003) is a flux limited sample of U /LIRGs from the 
IRAS All Sky Survey. Of the 629 extragalactic objects 
with 60 /urn flux greater than 5.24 Jy, 86 are interacting 
galaxies that are visually similar to our sample in the op- 
tical (i.e. close binaries). Of these, 32 have high-quality 
X-ray data from either Chandra or XMM-Newton that is 
sensitive to the presence of an AGN. Not accounting for 
the presence of undetected Compton-thick nuclei, only 
3% (1/32) of the RBGS sources with X-ray data show 
binary X-ray nuclei (NGC 6240; Komossa et al. 2003). 
This is consistent with the 0-8% (0-1 out of 12) we ob- 
serve in our modest SDSS sample. 

4.4. Frequency of Binary AGNs in SDSS Mergers 

From the very short snapshots of our study, we have 
found that coincidence of optically selected active nucleus 
with mildly obscured X-ray nucleus is relatively common 
(8/12). Given the faint detections, these snapshots are 
too short to place strong limits on the absence of AGN 
in the undetected galaxies, so it is difficult to comment 
on the frequency of binary active nuclei. However, we do 
detect a pair of X-ray nuclei in GZ 9, implying that this is 
uncommon unless the second nucleus is heavily obscured. 
In that instance, the most likely scenario would be that 
all nuclei are obscured. That is, either binary nuclei are 
uncommon, or merger nuclei in general have a high prob- 
ability of being heavily obscured. The latter possibility 
cannot be addressed by the current sample. To do bet- 
ter, we will need to increase the exposure times, expand 
our merger sample for better statistics, and include a 
sample of major mergers for which there are no optically 
detected nuclei. Another natural follow-up would be to 
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extend the study to a similarly select sample with a lower 
mass limit to examine the dependence of binary AGNs 
on the mass of the host galaxies. 
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TABLE 1 

The Sample and Results 
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Note. — Col.(l): Galaxy Zoo morgor identifier in this paper. * represents the nucleus in the pair with an optical spectrum from SDSS- Col. (2): SDSS spectroscopic object ID with, a prefix of 5877. Col.(3)-(4): Right 
ascension and declination in J2000- Col. (5): Redshift. Col. (6): Galactic column density in units of 10“*® cm - ^(Dickey ic Lockmau 1990). Col. (7): Projected optical separation of the nuclei in arcminutes (kpc). Col. (8) 
Optical spoctral type classification based on Figure 1 (A=AGN, L — LINER, C = AGN/Star-forrning composite, and S = ..tar-forming). Col. (9): log of the galaxy mars, derive 1 from the SDSS data (e.g., Schawinski et a|. 2010). 
Col.(10)-(11) : Number of counts in the soft (0.5-2 kcV) and hard (2—8 keV) bands, i-ejpectively. Co). (12): HardnotJ ratio. The error bars are propagated asaumiug Poisson statistics (Gehrels 1986). Col. (13): Photon index 
estimated from the observed hardness ratio assuming Galactic Njj. Col. (14): Estimated intrinsic column density for the source assuming a power law with V — 1.8 in units of 10^® on"*. Col. (15): Estimated 2-10 keV flux 
by asaumiug a rodshifted power law with utimated T from Column 13 normalised by the detectod count rate. The flux is given in units of 10“ 14 erg s“ 1 cin“ 2 . Col. (10): Estimated 2-10 keV luminosity based on flux in 
Column 1ft. The luminosity Is given in units of 10*® erg s — * . 
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Fig. 1. — Emission line diagnostic diagram used for the selection of the Chandra Galaxy Zoo sample. The gTey scale represents the 
complete Galaxy Zoo sample and the green points are the Galaxy Zoo AGN selected using narrow line diagnostics from Schawinski et al. 
(2010). The dashed curve shows the empirical separation between purely star-forming galaxies and the composite region of the diagram as 
determined by Kauffman et al. (2003, Ka03). The solid curve is the theoretical extreme starburst line of Kewley et al. (2001, keOl) beyond 
which the dominant source of ionization must be due to something other than star formation. The straight line demarcates the empirical 
AGN-LINER separation in Schawinski et al. (2007, S07). The yellow points are the Galaxy Zoo mergers that meet the criteria of mergers 
having mass ratios of at least 3:1 and having one of the progenitors with a mass > 10 11 ATg with significant emission line detections; 
most are LINERs. The blue points are the secondary nuclei with SDSS optical spectra. The Chandra observed nuclei (both primary and 
secondary) are circled in red with their GZ identifier from Table 1 labeled. Our Chandra sample covers the full range of [O III]/H/8 and 
[N II]/Hq emission line ratios for composite and LINER objects. They are also representative of the merger distribution in the composite 
and LINER areas of this Baldwin-Phillips-Terlevich (BPT) diagram (Baldwin et al. 1981). 



Fig. 2, — Composite gri images of the Chandra-observed sample from SDSS Data Release 7. Each frame is labeled with the Galaxy Zoo 
identification listed in Table 1 and the horizontal bar represents angular distance of 20 arcseconds. 





Fig. 3. — False color X-ray images of two interesting objects in our sample, GZ 9 (left) and GZ 10 (right). The red represent the 0.5-2 keV, 
green the 2-6 keV, and blue the 6r-8 keV emission. The contours are from SDSS i band images. The raw X-ray images were smoothed with 
a Of'5 Gaussian, the width of the nominal point spread function of Chandra. The binary nuclei in GZ 9 are both detected in the X-ray, 
though the southern nucleus is dominated by soft X-ray emission. In GZ 10, the X-ray emission shows east-west extension which may be 
due to star formation in addition to an obscured AGN. The flatness of the X-ray spectrum implies a column density ~5x 10 21 cm -2 if 
we assume a canonical power-law photon index of F ~ 1 .8. 
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Fig. 4 - Comparison of star formation rates derived from the 2-10 keV luminosity (Ranalli et al. 2003) and SDSS u band luminosity 
density (Hopkins et al. 2003). Only the detected nuclei with hardness ratios, and thus estimated X-ray luminosities, are plotted. Errors 
are 1 — <r. The arrows indicate poorly constrained negative error bar for three nuclei whose HR lower limits approach -1, where the Nh 
and F values become degenerate for any value of HR. The solid line is the line of equal star formation rates. Each of the detected nuclei 
is labeled corresponding to the identification in Table 1. All of the detected nuclei have X-ray luminosity above that expected from star 
formation. 
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Fig. 5. — Cumulative rest-frame photon distribution in the hard band (2-8 keV) for the 12 detected nuclei. The black soiid histogram 
represents the combined detected counts from all 12 detected sources in our sample. The blue dotted curve represents the expected 
cumulative distribution assuming an unobscured AGN model where a single unabsorbed power law represents the source of the AGN 
emission. The red dashed curve shows the distribution for a Compton-thick AGN model where the emission is represented by a power law 
plus a 6.4 keV iron emission line with an equivalent width of 1 keV. White there is a clear difference between the two modeled distributions, 
when compared with the observed distribution, neither can be shown to be the preferred model by a two-tailed K-S test. Therefore, we 
cannot rule out contribution from Compton-thick AGNs at present. 




